Two-dimensional fermionic superfluids, pairing instability and vortex liquids in the 

absence of time reversal symmetry 
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We consider a generic two-dimensional system ol lermionic particles with attractive interactions 
and no disorder. II time-reversal symmetry is absent, it is possible to obtain incompressible insulating 
states in addition to the superfluid at zero temperature. The superfluid-insulator phase transition is 
found to be second order in type-II systems using a perturbative analysis of Cooper pairing instability 
in quantum Hall states of unpaired fermions. We obtain the pairing phase diagram as a function 
of chemical potential (density) and temperature. However, a more careful analysis presented here 
reveals that the pairing quantum phase transition is always preempted by another transition into a 
strongly correlated normal state which retains Cooper pairing and cannot be smoothly connected 
to the quantum Hall state of unpaired fermions. Such a normal phase can be qualitatively viewed 
as a liquid of vortices, although it may acquire conventional broken symmetries. Even if it did not 
survive at finite temperatures its influence would be felt through strong quantum fluctuations below 
a crossover temperature scale. These conclusions directly apply to fermionic ultra-cold atom systems 
near unitarity, but are likely relevant for the properties of other strongly correlated superfluids as 
well, including high temperature superconductors. 



Many years of efforts to understand the strongly cor- 
related phases of cuprate superconductors and other un- 
conventional materials have resulted in a variety of new 
conceptual ideas. One broadly accepted idea is that the 
destruction of two-dimensional superconductivity in un- 
derdoped cuprates involves strong fluctuations of the or- 
der parameter instead of Cooper pair dissociation 1 . The 
remarkable theoretical insight gained from duality trans- 
formations in several systems 2 -*^ stimulates viewing this 
phase transition as a proliferation of vortices. The nor- 
mal state may be envisioned as a liquid of mobile vor- 
tex loops. Indeed, there are several experimental obser- 
vations that support this view 4 -! 5 -, and several theories 
that exploit ft 6 -' 7 - to successfully explain certain uncon- 
ventional properties of cuprates^ 9 - ' 10 ' 11 ' 12 . 

Nevertheless, the notion of a vortex liquid is not a triv- 
ial one. It is not even clear yet whether any correlated 
state of any known material truly deserves to be called 
a vortex liquid. In particular, the vortex liquid regime 
found in the pseudogap region in cuprates seems to be 
separated from the high temperature normal state only 
by a crossover. A quantum vortex liquid might be a 
genuine phase at zero temperature which influences the 
finite temperature dynamic s 13 ' 14 . A few theoretical mod- 
els can reliably describe strongly correlated phases that 
can be characterized as vortex liquids, and they provide 
the best route to precisely define and characterize these 
phasesii 1 ^. 

A new platform for addressing these interesting chal- 
lenges in many-body physics has been recently intro- 
duced in atomic physics 17 . The ability to condense 
ultra-cold fermionic atoms into a superfluid with control- 
lable strength of interactions and study them as virtu- 
ally ideal and clean systems has already mobilized an ef- 
fort to simulate condensed matter systems. Static three- 
dimensional cold-atom systems near a broad Feshbach 
resonance have striking universal propertie s 1 i 18 ' 19 ' 20 



which make them amenable to field theoretical analysis. 
This universality is controlled by a quantum critical point 
at the zero density Feshbach resonance^ 9 -, which defines 
the unitarity limit. The superfluid in this regime interpo- 
lates between the Bardeen-Cooper-Schrieffer (BCS) state 
of Cooper pairs and Bose-Einstein condensate (BEC) 
of diatomic molecules, and insulating states in optical 
lattices interpolate between the two analogous limits of 
band and Mott insulators respectively^!. The cold-atom 
superfluids can be destabilized by fast rotation in the 
same way the electronic superconductor can be desta- 
bilized by magnetic fiel d 22 ' 23 ' 24 . Experimental progress 
toward observing this effect in cold atoms is still limit- 
ted, but the overall pace of development in the field gives 
many reasons for optimis m 25 ' 26 ' 27 . 

In this paper we analyze a simple two-dimensional 
model of interacting fermionic particles which is a signifi- 
cant simplification of electronic systems such as cuprates 
in magnetic field, but can be routinely and accurately 
realized using rotating ultra-cold atoms. The model 
contains paired superfluids at low temperatures and a 
normal phase of uncorrelated fermions at high temper- 
atures. We concentrate on type-II superfluids which 
host an Abrikosov vortex lattice due to the absence of 
time-reversal symmetry. The uncorrelated normal phase 
is a thermally excited quantum Hall state of unpaired 
fermions. A characteristic phase diagram of pairing in- 
stability in this normal phase is shown in FigflJ One of 
our main results is that the normal phase can stretch all 
the way to zero temperature at finite particle densities 
for certain strengths of interactions between fermions. 

The pairing instability depicted in Fig[T]was examined 
in the past by several authors; for a review see 2 ^ and 
references therein. While three-dimensional systems at- 
tracted most attention in the literature, the focus here 
is on two-dimensions precisely because of the previously 
missed zero-temperature normal states (they are not ex- 
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FIG. 1: Critical temperature of the superfluid-normal pairing 
transition as a function of chemical potential. Interactions be- 
tween particles are characterized by the two-body scattering 
length a and become weaker going from (a) to (d): v oc — aT 1 
(defined in the text and Fig(3} is 0.3 in (a), 0.4 in (b), 0.5 in 
(c), 0.6 in (d). The vertically stretching lines are normal state 
constant density contours, with increment 0.25B/(2n) where 
B = 2mu is the uniform flux density. 



pected to occur in three dimensions). In this paper we 
provide a fresh approach to issues such as superconduc- 
tivity in high magnetic fields and keep the formalism sim- 
ple and clear in order to focus on the most fundamental 
questions. This simplicity will allow us to reach some fun- 
damental conclusions which were obscured in the past by 
the inability to go beyond the mean-field and other ap- 
proximations. We will not consider Zeeman effect in this 
paper; Zeeman pair-breaking effects can be completely 
eliminated in cold-atom systems, while in electronic sys- 
tems they are certainly important although the super- 
conducting state may be able to survive them even in 
high magnetic fields (this analysis will be published sep- 
arately). Another popular approach in the literature has 
been the Landau- Ginzburg theory, but for our purposes 
it is important to take into account the internal structure 
of the bosonic degrees of freedom responsible for superflu- 
idity and capture phenomena characteristic for fermionic 
superfluids that are absent in pure bosonic systems and 
Landau-Ginzburg theories. 

The possibility of finding a zero-temperature normal 
phase of unpaired fermions, which is actually a quan- 
tum Hall insulator, provides an opportunity to estab- 
lish with physical rigor the existence of another non- 
superfluid phase at zero temperature. We show that the 
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FIG. 2: (color online) A schematic phase diagram at zero 
temperature. The depicted phases are fermionic quantum 
Hall insulator (QHI), vortex liquid (VL) and superfluid with 
an Abrikosov vortex lattice (SF). The vertical axis is a mea- 
sure of the attractive interaction strength between fermions. 
The horizontal axis in is fermion chemical potential expressed 
in reference to a Landau level (w c is cyclotron frequency; the 
left and right edge of the diagram correspond to two neigh- 
boring Landau levels). In the weak pairing limit the super- 
fluid can be destabilized between Landau levels, forming a 
vortex liquid across a first order phase transition. For very 
weak pairing interactions a quantum Hall insulator of un- 
paired fermions might become stable across a transition such 
as the one discussed in Ref.™. The particle densities in VL 
and QHI are expected to be fixed and commensurate with 
respect to magnetic area. In this paper we find integer QHI 
states shown in Fig[TJand the corresponding VL states around 
them. Other insulating states such as fractional bosonic and 
fermionic quantum Hall states can be expected in various cir- 
cumstances. 



zero-temperature superfluid at a finite density must be 
destroyed by the quantum fluctuations of the order para- 
meter before the depairing instability occurs. This vor- 
tex lattice melting quantum phase transition is expected 
to be first orde r 29 ! 30 ' 31 . The resulting normal phase is 
strongly correlated as it retains Cooper pairs, and it can- 
not be smoothly connected with the quantum Hall state 
of unpaired fermions. Quantum phase transitions be- 
tween the two kinds of normal phases have been stud- 
ied recently^. We will refer to the new normal phase 
as vortex liquid, based on its sharp distinction from the 
fermionic quantum Hall states. A schematic phase dia- 
gram is shown in FigEJ 

The concrete properties of a vortex liquid and its sta- 
bility at finite temperatures are not addressed in this 
paper. The main reason is that these properties are 
not universal in generic realistic systems. To make this 
point we use a renormalization group analysis and show 
that all possible interactions are relevant at the Gaus- 
sian fixed point in two dimensions, regardless of their 
spatial dependence or the number of colliding particles. 
The existence of infinitely many relevant directions in the 
parameter space suggests the possible existence of a large 
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number of stable interacting fixed points which determine 
the properties of phases. Which of these phases is real- 
ized in a particular realistic system is highly sensitive to 
the microscopic details of interactions between particles. 
One prominent candidate in the case of short-range in- 
teractions is a quantum Hall liquid of incoherent Cooper 
pairs 24 . Other possibilities include a conventional bro- 
ken symmetry such as a density wave, in which case its 
algebraically correlated analogue could exist as a stable 
phase at finite temperatures. 

Since the ultra-cold fermionic atoms near unitarity are 
the closest realization of the model considered here, we 
will devote extra attention to their specifics and express 
the results in scales natural to cold atoms. However, the 
analysis is general and possibly relevant to other strongly 
correlated systems such as cuprates. 

This paper is organized into several sections. We begin 
by defining the model of interest in section fX] and set- 
ting up the formalism for perturbation theory. In section 
fBl we analyze the phase diagram of the Cooper pairing 
transition in the perturbation theory. The interpretation 
and limitations of the perturbative results are outlined in 
section [C] In section [D] we establish that vortex liquid 
phases, missed by the naive perturbation theory, exist at 
zero temperature. At the end, we apply a renormaliza- 
tion group analysis in section [^] to appreciate the non- 
universal properties of the vortex liquid states and the 
richness of its possible physical realizations. We sum- 
marize the results and discuss a possible experimental 
approach in the concluding section [Fj 



The model 



when interactions are strong. The physical case of inter- 
est is N = 1. The attractive interactions between parti- 
cles are decoupled in the particle-particle channel by the 
Hubbard-Stratonovich field $. These interactions lead to 
processes in which two colliding particles are annihilated 
and recreated at a different location and different time, 
so that we can view $ as a mediating Cooper pair field. 
No particular assumptions are made about the nature or 
range of interactions at this point, so we leave the bosonic 
kernel Ho unspecified until the next section. 

It is worthwhile noting that in cold-atom systems the 
gauge field A is rigid and cannot be expelled from a su- 
perfluid because it captures the inertial forces in the ro- 
tating frame of reference. As a consequence, cold atom 
superfluids are automatically type-II. In contrast, A is 
dynamical in electronic systems and the type of a su- 
perconductor is determined by various parameters not 
contained in this model, such as disorder. 

The action is formally a generalization of the "two- 
channel" model to the Sp(N) symmetry group. The pur- 
pose of this generalization is to systematically organize 
a perturbation theory with l/N as a small parameter 
even when interactions are not weak, as will be the case 
throughout the paper. We first note that the full Cooper 
pair propagator is proportional to 1/N. This follows from 
the Dyson equation whose solution can be represented 
using Fcynman diagrams: 
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Our starting point is a generic model of fermionic par- 
ticles with attractive interactions in two dimensions. We 
couple the particles to a U(l) gauge field which imple- 
ments a uniform flux density B: 



V x A = zB = zmu> c — zlmto 



(1) 



In electronic materials flux is created by the external 
magnetic field B, while in ultra-cold atom systems it orig- 
inates in rotation at angular velocity u> (the relationship 
above follows from dynamics in the rotating frame of ref- 
erence). We will not consider the fluctuations of A, so 
that the model will naturally describe neutral particles 
(although Coulomb forces can be implemented as direct 
interactions). The imaginary-time action is: 
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The fermionic (Grassmann) fields if}i a carry spin a G {] 
, 1} and a flavor index i = 1 . . . N which will be used 
to systematically generate perturbative expansions even 



A double line represents IL7 1 , while a single thick line 
is the full fermion propagator (note that the summation 
over fermion flavors in the bubble diagram contributes a 
factor of N) . The full vertex represented by a shaded dot, 
and the full fermion propagator are obtained by dressing 
the bare vertex and propagator with full Cooper pair lines 
from the left-hand-side, and this generates contributions 
at higher orders in 1/N. In general, any physical observ- 
able will formally be expressed as an expansion in powers 
of 1/N, justifying perturbation theory for any interaction 
strength in the limit of large N. This is a semi-classical 
expansion because for N = oo all pairing fluctuations 
are suppressed and saddle-point approximations become 
exact. 

The uncorrelated Gaussian state to which the pair- 
ing interactions are added as perturbations is an integer 
quantum Hall state of fermions. This state is hard to 
imagine as a ground state of any interacting system be- 
cause of its macroscopically Landau-degenerate spectra. 
Nevertheless, it is a good starting point at sufficiently 
high temperatures, and its usefulness at zero temperature 
will become apparent in section [Dj 

The bare fermionic states are Landau orbitals with en- 
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TABLE I: Bare elements of the perturbation theory in Landau 
representation. The vertex function V is given by (O, and 
ui m is Matsubara frequency. 



ergies e n which we will represent in the Landau gauge: 

A = -Byx (4) 

( 1 

e« = w c I n + - 

The quantum numbers are Landau level index n € 
{0, 1,2...} and momentum k in x-direction. The bare 
fermion wavefunctions an d the Cooper pair wave- 
functions & n ,p that will be used to expand the order 
parameter are: 
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(H n are Hermite polynomials). Note that $„ p carries 
twice the "charge" of fermions ip n k- 

The elements of the perturbation theory in Landau 
representation are summarized in table [T] The dimen- 
sionless vertex function in Landau representation is read- 
ily derived from (|5|): 
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Note that the vertex does not depend on the momentum 
transferred to the Cooper pair, but only on the difference 
of incoming fermion momenta. 



B. Pairing instability 

Consider bosonic two-particle excitations in the nor- 
mal phase. Their dynamics is captured by the Cooper 



pair propagator D = (NIl)^ 1 . In general, we need to 
view this propagator as an operator, or a matrix whose 
rows and columns are indexed by the quantum num- 
bers of Cooper pairs in some representation. We do not 
know a priori which representation diagonalizes II, but 
the Landau representation is a convenient one to work 
with. Since pairing instability occurs at zero frequency, 
the easiest way to detect it is by the onset of negative 
eigenvalues of II at zero Matsubara frequency. This fol- 
lows from the fact that upon integrating out the fermion 
fields in © the effective Cooper pair action has the 
quadratic term $tll<j> which needs to be positive for any 
vector $ if the perturbation theory about the normal 
state ($ = 0) were to be stable. 

Before embarking on the calculation of II we must note 
that it is not gauge invariant. Gauge transformations 
in position representation that leave the action invariant 
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Gauge transformations are manifestly unitary transfor- 
mations of II and do not affect its spectrum of eigenval- 
ues. Therefore, the spectrum which contains all informa- 
tion about bosonic dynamics is gauge invariant and the 
calculations can safely proceed in any gauge choice. Note 
that if the gauge field were dynamic one would not be 
able to apply the standard perturbation theory in pow- 
ers of charge because the gauge transformations would 
appear non-perturbative. That problem is absent in our 
approach. 

The Cooper pair kernel II can be viewed as a sum of 
two parts according to §5§ . The first part is the bare ker- 
nel IIo which contains information about the microscopic 
details of the attractive interactions between fermions. 
The second part is the fermion bubble diagram which 
can be derived as an expansion in 1/N: 
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Note that this expression for IT' written in the Landau 
representation is diagonal in p but not in the bosonic 
Landau level indices n. The lowest order term, explic- 
itly written here with the internal Matsubara frequency 
summed up, contains only the bare fermion propagators 
and no vertex corrections. 

As we will discuss in section [j3 different microscopic 
details of attractive interactions between fermions can in 
principle lead to different properties of the system. With- 
out being able to classify all possible cases we still have 
to choose a particular form of interactions in order to 
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carry out calculations. The simplest choice is to model 
short-range interactions by a zero-range potential. The 
Hubbard-Stratonovich decoupling of such attractive in- 
teractions produces ITo = const.: 



exp | J drd 2 rV4 T (r )V4 X 0) (r)^ n (r) 
X^expj-^ drd 2 r 



(9) 
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In three dimensions and zero flux density the constant 
ITo = V^ 1 is related to the scattering length of two body 
collisions in vacuum: ITo = mis/An, where a = —\jv is 
the scattering length. In the context of cold atom exper- 
iments, v is the detuning from the Feshbach resonance. 
For v — the scattering length diverges and the system 
acquires universal properties, which is the so called uni- 
tarity regime where microscopic details of interactions do 
not matter and the simple interaction potential above is 
sufficient for making various quantitative predictions. In 
two dimensions we must redefine the detuning in order 
for IIo to have proper engineering dimensions: we re- 
place v by v/a z , where a z is a confinement length-scale 
or "width" of the potential well in z-direction which ef- 
fectively reduces the system to two dimensions. 

While establishing these connections with cold-atom 
systems is useful, one should be aware of the conceptual 
problems in the presence of finite flux density. Strictly 
speaking, scattering length a is a quantity that describes 
interactions between extended free particles, but in our 
case particles are localized in Landau orbitals. Only for 
a <C R c , the cyclotron radius at Fermi energy, a retains 
its physical meaning. Therefore, it is most appropriate 
to regard v = —l/a just as a measure of the strength 
of interactions, positive values being the BCS limit and 
negative values being the BEC limit. Interpreting v = 
as a Feshbach resonance in cold-atom systems and as- 
cribing universal properties to it is problematic, and the 
discussion in section fEl will shed more light on this. 

In fact, the relationship between IIo an d scattering 
length is established only after the regularization of the 
infra-red and ultra-violet behavior of II. In three dimen- 
sions and without the external flux the II given by the 
expression analogous to ([8]) is known to be ultra-violet 
divergent. Once this field-theoretical artifact is removed 
by regularization, the leftover contribution to II is pre- 
cisely mv/4ir— . In two dimensions II becomes infra-red 
divergent as well if the bare fermion spectrum is gapless. 
In our case the spectrum is gapped due to Landau level 
quantization, but more generally any two-dimensional 
system should be viewed as a three-dimensional system 
constrained by a potential well in which the lowest energy 
band is raised above the bottom of the well. Therefore, 
infra-red divergence is not a real concern. The ultra- 
violet behavior is, however, very different from the zero 




FIG. 3: (color online) Critical temperature for the onset of 
Cooper pairing as a function of chemical potential for several 
values of detuning v — v j '(4mtja z ): the outermost curve is 
BEC limit v = -0.9, the innermost is BCS limit v = +0.9, 
unitarity v = is dashed (Au = 0.1). The vertically stretch- 
ing lines are normal state constant density contours, with in- 
crement 0.25B/(2tt). 



flux case, because the flux modifies the spectrum at ar- 
bitrarily high energies. On physical grounds we perform 
regularization by subtracting from ([5]) the same expres- 
sion II' (0) evaluated at zero Matsubara frequency and 
momentum transfer, as well as zero temperature and 
chemical potential, but finite flux density. This takes 
care of any ultra-violet divergences and is compensated 
by adding mv/Anciz which in the limit of small scattering 
length is precisely given by n'(0). The final expression 
for the Cooper pair kernel that we shall use is: 



n 
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n' - n'(o) + o 



(10) 



Figure [3] shows pairing second-order phase transitions 
in the limit N — > oo obtained from the onset of negative 
eigenvalues of (fl~0"|) . The numerical calculation of (JS]) 
involves the first 500 fermionic Landau levels below the 
cut-off and the lowest bosonic Landau level, although the 
phase diagram remains the same if higher bosonic Lan- 
dau levels n, n' are included. The deeper one goes into the 
BEC regime, the faster the growth of critical temperature 
with chemical potential (density). In the BCS regime, 
however, the paired phase breaks up into dome-shaped 
islands at low densities, obtained when chemical poten- 
tial gets close to a fermionic Landau level. Once a paired 
state becomes a connected region at sufficiently large 
densities, there are trailing islands of unpaired states sit- 
ting between the Landau levels. These details are shown 
with more clarity in Fig[TJ Related phenomena have been 
considered in literature, but the zero temperature normal 
states were not anticipated^. 
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Integer quantum Hall states are obtained by putting 
the chemical potential \i between Landau levels. Pair- 
ing interactions compete with the cyclotron gap lo c = 2uj 
and must be strong enough in comparison to the small- 
est single-particle gap min{|/i — u> c (n + h)\, n 6 Z} in 
order to give rise to Cooper pairing. This means that 
if pairing is not too strong, insulating states can exist 
at zero temperature with fully populated Landau levels 
below the chemical potential. The analysis based on pair- 
ing instability in the perturbation theory is, however, not 
complete. In the next two sections we discuss what has 
been missed. 



C. Limitations of the perturbation theory 

The perturbation theory seemingly discovers Cooper 
pairing instability for any finite N when the chemical 
potential is brought sufficiently close to a Landau level. 
To cure this instability one usually redefines the ground- 
state by introducing an order parameter and then re- 
builds the perturbation theory about the new ground- 
state. The order parameter which describes a superfluid 
with broken U(l) symmetry and an Abrikosov lattice of 
vortices can be written as a superposition of Landau level 
wavef unctions in ([5]): 

$ o(r) = E / Tr^WO • (11) 

Hence, the order parameter is specified by multiple am- 
plitudes 4> n ,p- If we substitute <&(r) = &o(r) + 6$(r) 
in and integrate out the fluctuating fields ip and <5$ 
we can obtain the perturbative I/TV expansion for free 
energy density of the ansatz state characterized by <I>o : 

^1 = ^ + ii,,«!>;/.f;; + u m H*H** k < + om . 

The quantum numbers (n,p) are represented by indices 
i,j... for brevity, and all couplings IL, U . . . are func- 
tions of fx,T,u> . . . which have tensor structure and con- 
tributions to all powers of 1/N. Note that the expansion 
is analytic at 3>o = because the unperturbed state at 
$o = is a fully gapped incompressible quantum Hall 
state of unpaired fcrmions. In the normal state T is min- 
imized at $o — 0, while in the superfluid the minimum 
is at $ 7^ 0. 

At the lowest order or the mean-field level we found 
both superfluid and normal states in the previous sec- 
tion. If we include higher order corrections the order 
parameter in a superfluid state will be reduced and the 
normal areas in the phase diagram will grow. However, 
superfluid states will survive at all orders of perturba- 
tion theory for any finite N if the chemical potential is 
sufficiently close to a Landau level. This is ultimately a 
consequence of the cyclotron gap in the fermion spectrum 
which prevents the occurrence of non-analytic infra-red 
features and yields finite corrections to all couplings in 



(fT2t at all orders of 1/N. Alternatively, any point in the 
phase diagram which is a superfluid in the mean-field ap- 
proximation will remain superfluid for sufficiently large 
N despite fluctuations. 

Throughout this paper we exploit the computational 
and analytical convenience of the physical picture embod- 
ied in (|12p . However, we must be aware of its limitations. 
This picture naively suggests two kinds of broken sym- 
metries in the superfluid phase, the "off-diagonal" long- 
range order and space-group symmetry breaking due to 
the vortex lattice. Even at finite temperatures it sug- 
gests the existence of these broken symmetries. On the 
other hand, it is well known that fluctuations in two- 
dimensional systems can restore continuous symmetries 
in the equilibrium states. In particular, the superfluid 
phase at zero temperature has a long-range ordered vor- 
tex lattice, but shear fluctuations of vortex positions re- 
duce the "off-diagonal" order to algebraic correlations 
and hence restore the U(l) symmetr y 32 i 33 ' 34 . At finite 
temperatures even the vortex lattice must be melted ac- 
cording to Mermin- Wagner theorem, but algebraic cor- 
relations of vortex positions are allowed and Kosterlitz- 
Thouless transitions can separate the quasi long-range 
ordered vortex "lattice" phases from disordered phases 
with only short-range correlations. 

The naive approach behind the Eqfl2] assumes that 
the only two kinds of states one needs to worry about 
are the normal state with short-range correlations of the 
order parameter, and the superfluid with long-range cor- 
relations. It does not explore the free energy density T' 
of states with algebraic correlations of the order para- 
meter (s). It turns out that in the circumstances men- 
tioned above T' < J-(Qq) for any attempted order para- 
meter $o inside the naively predicted superfluid regions. 
Another possibility is that T' < ^(0) inside some parts 
of the naively predicted normal states. 

Calculating T' is highly non-trivial in the perturbation 
theory. Therefore, we are better off learning what we can 
from JF($o) in 02). The difference A = JF($ ) - F de- 
fines an energy scale associated with the error made by 
(fT^)l . We will be able to make conclusions only about 
phenomena characterized by energies larger than A, or 
length-scales smaller than S cx A -1 . This will be suffi- 
cient for establishing the existence of zero temperature 
vortex liquid phases in the next section. On general 
grounds we can expect that the error A is small because 
even in the absence of a true long range order the scale 
given by $o plays a role in various properties of the sys- 
tem. 



D. Vortex liquid 

The quantum Hall insulator of unpaired fermions is a 
very peculiar unperturbed state for building the pertur- 
bation theory. A very important property of the exact 
Cooper pair propagator D — (A^II) -1 in this state is 
that its spectrum is macroscopically degenerate and cor- 
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FIG. 4: Feynman diagram labeling examples. 



responds to bosonic Landau level degeneracy. A simple 
argument is that in the absence of any emerging param- 
eters that could characterize the equilibrium state the 
propagator D obtains its matrix structure from the only 
single-particle operator available in the theory ((2]), the 
canonical momentum P = — iV — A. This operator has 
a degenerate spectrum that gives rise to Landau level 
degeneracy, hence the operator D(P) must also have a 
degenerate spectrum, assuming of course that the func- 
tional dependence on P is analytic. A detailed demon- 
stration of this fact follows from perturbation theory to 
all orders, and agrees with a general expectation that 
the perturbation theory describes a state smoothly con- 
nected to the state of non-interacting particles in the ab- 
sence of instabilities. It is important to note here that 
Cooper pairs are "charged" with respect to the gauge 
field; "charge-neutral" particle-hole excitations are ex- 
empt from the Landau-level quantization. 

The perturbative argument goes like this. We can es- 
tablish that ijHJ) and hence (TIT)]) are independent of p to 
all orders of 1/N. First, the written lowest order term in 
the 1 /N expansion involves fermion energies that do not 
depend on ±k + p/2, and vertices (0 that do not de- 
pend on p. This term defines the bare boson propagator 
Dq oc N , using which we generate higher order correc- 
tions of ©. Let us apply the following labeling rules 
in all Feynman diagrams that contribute to (|10[) : each 
fermion propagator shall carry momentum k+p/2 in the 
arrow direction, while each boson propagator shall carry 
q + p (see FigH]). The transfer of p is automatically con- 
served. Each added vertex takes momenta k\ +p/2 and 
k'2 +p/2 at its fermionic terminals, but there is no depen- 
dence on p since only the difference of momenta at the 
fermion terminals matters in ^ . The added bare boson 
lines take p contributions, but according to ([5]) they do 
not depend on the transferred momentum. Therefore, 
order by order, no dependence on p is introduced in any 
Feynman diagram. 

Changing the nature of interactions does not modify 
this conclusion because Ho must be gauge invariant in 
the presence of the external uniform flux density. Even 
interactions in the particle-hole channel cannot lift the 
degeneracy of bosonic modes (an equivalent "degeneracy 
theorem" can be proved to all perturbative orders in a 
Landau-Ginzburg theory of interacting particles in mag- 
netic field). Therefore, the gapped two-particle bosonic 
modes are macroscopically degenerate in the trivial insu- 
lating state regardless of the type of interactions between 
particles. Naturally, it is hard to believe that any quan- 



tum fluctuations would not lift this degeneracy. At no 
surprise we shall demonstrate that the degeneracy must 
be lifted at least close to the superfluid transition, but 
the main point for now is that it takes spontaneous re- 
organization of the ground-state to lift the degeneracy 
(which in some sense is a trivial but very important ob- 
servation). Even without knowing whether the obtained 
normal T — state has some conventional order or topo- 
logical order, it is certain that it cannot be smoothly 
connected with the ordinary insulating state. 

First, we need to understand the properties of the zero 
temperature superfluid phase. Let us start from a quan- 
tum Hall insulator at zero temperature. When the energy 
of the lowest degenerate bosonic states in the insulator 
becomes negative in response to the chemical potential 
approaching a Landau level, then Cooper pairs condense 
into those states. The resulting order that minimizes free 
energy density is automatically an Abrikosov vortex lat- 
tice residing initially just in the lowest bosonic Landau 
level. The free energy cost of rearranging vortices with- 
out changing the vortex density is exceptionally small 
near the phase transition. By assumption we are dealing 
with type-II superfluids so that the pairing phase transi- 
tion is second order. But, even more striking is the effect 
of bosonic degeneracy in the trivial insulating phase. 

Consider the free energy density as a function of or- 
der parameter in (TT2|) . We argued that the quadratic 
coupling fl is macroscopically degenerate. Therefore, dif- 
ferent vortex lattice arrangements, specified by different 
order parameters will differ in free energy only as 
0(<i>o) close to the phase transition. We conclude that 
the potential energy density of a vortex lattice depends 
on the order parameter near the transition as: 

Spot ~ $0 ■ (13) 

Quantum fluctuations of the order parameter dynam- 
ically distort the vortex lattice. Their microscopic ori- 
gin are the short-range repulsive interactions between 
Cooper pairs which are caused at least by the Pauli exclu- 
sion of their fermionic constituents. Interactions between 
bare fermions also affect the quantum fluctuations in a 
manner that can be very sensitive to microscopic details 
of the interaction potential. In fact, all possible two- 
body and multi-particle interactions are relevant at the 
Gaussian fixed point in the renormalization group sense, 
regardless of their spatial dependence (see next section) . 
Having infinitely many relevant directions in the para- 
meter space suggests that there could be many interact- 
ing stable fixed points that describe numerous possible 
correlated phases. Systematically analyzing all possibil- 
ities seems hopeless. Therefore, we shall not attempt 
to identify the properties of strongly correlated normal 
phases that can be obtained after the superfluid is de- 
stroyed by quantum fluctuations. Instead, we can at least 
prove that some strongly correlated normal phase must 
inherit the superfluid, rather than the trivial unpaired 
quantum Hall insulator. This is already one statement 
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about the renormalization group flow of parameters near 
the superfluid quantum phase transition. 

In order to summarize the effects of quantum fluctua- 
tions we can view them as quantum motion of vortices. 
In this picture, vortex lattice melting destroys the super- 
fluid. It will suffice to estimate the effective vortex mass 
M v , specifically its dependence on the order parameter 
$o- We make the estimate using the Heissenberg uncer- 
tainty relation. Vortices are localized in a vortex lattice 
at the scales of magnetic length l m = 1/ ^Jmuj c which 
measures the average distance between vortices and cor- 
responds to a momentum scale p m ~ ^/rmu c . If vortices 
were allowed to move freely they could convert the poten- 
tial energy of localization in the lattice to kinetic energy 
of the order of p^ n /2M v per vortex. Since vortex den- 
sity is topologically fixed by the constant rotation rate 
or magnetic field, we find: 

£ kln ~ M- 1 ^,,) • (14) 

In typical fermionic superfluids one is usually mostly 
concerned with the quasiparticle contribution to vortex 
mass. In s-wave superfluids vortex cores can localize 
quasiparticles, so that M v is at least of the order of the 
total mass of trapped core fermions. An upper bound for 
M v can be found by viewing a vortex core as a potential 
well of depth | $o I in the units of energy, and radius com- 
parable to the Cooper pair coherence length £ ~ |$o| _1 - 
Assuming roughly a constant density of core states ap- 
propriate for two-dimensional non-relativistic quasipar- 
ticles, the total number of core states is proportional to 
| $o | , so that M v < | $o | • This estimate is best suited for 
the strong BCS limit 3 ^. As the strength of interactions is 
increased, the coherence length decreases, quasiparticles 
are gradually being expelled from vortex cores and M v is 
being reduced. There is also a contribution of extended 
quasiparticles to vortex mass. In s-wave superfluids this 
contribution is thermally activated, but in d-wave super- 
fluid it gives rise to a vortex mass of the order of a fermion 
mass 36 . 

Close to the superfluid transition it turns out that the 
"hydrodynamic" vortex mass due to the phase fluctua- 
tions of the order parameter dominates in M v (for any fi- 
nite N). This contribution can be generally estimated as 
M v ~ e v /c 2 , where e v is energy of a static vortex and c is 
the speed of sound, or the superfluid critical velocity 36,38 . 
For a neutral superfluid with stiffness p s in the units of 
energy, c ~ \J p s /m while e v oc p s log(/ m /£) diverges loga- 
rithmically with the separation l m between vortices. The 
logarithmic behavior comes from the contribution of the 
circulating supercurrent phase gradients that extend be- 
yond the vortex core of radius £ freely up to the distance 
to the neighboring vortices (we assumed that the London 
penetration length is larger than the magnetic length l m ). 
Using p s ~ I'&ol 2 we find that the hydrodynamic contri- 
bution to the vortex mass behaves as M v ~ log | $o I near 
the phase transition for a fixed flux density. In charged 
superfluids the logarithmic factor is screened-ou t 36 ' 39 , so 
that the vortex mass approaches a constant as |$o| — » 0. 



Therefore, the vortex kinetic energy density is: 

E _ f llogl^oll" 1 > neutral superfluid 1 ^ 
111 \ const. , superconductor J 

Regardless of the concrete microscopic details and fac- 
tors that enter (fl3|) and (fT"5|) . close enough to the pair- 
ing transition the kinetic energy wins and vortex lattice 
must melt. The quantum melting occurs at a finite |$o| 
below which E^in > E pot , before the depairing transi- 
tion is reached. Both the kinetic and potential energy 
of vortices can be calculated in the perturbation theory, 
and the melting transition is expected to be first order 
according to the Landau-Ginzburg theory of phase tran- 
sitions. The resulting normal state can be regarded a 
strongly correlated vortex liquid in which Cooper pairs 
remain the most important low energy degrees of free- 
dom. The energy released into the motion of vortices in 
the vortex liquid is not sufficient to overcome the pairing 
gap |$q| m the fermion spectrum. 

By comparing the spectrum of bosonic excitations in 
this vortex liquid and the trivial quantum Hall insulator 
one easily concludes that the two normal phases cannot 
be smoothly connected. The bosonic modes in the super- 
fluid state close to the pairing transition are nearly degen- 
erate, so that many different vortex arrangements have 
very similar free energy densities. In the trivial quantum 
Hall insulator all of these modes become gapped and col- 
lapse to the same bosonic Landau level, leaving the occu- 
pied fermionic Landau levels to define the ground-state. 
On the other hand, vortex liquid is obtained by mixing 
the nearly-degenerate bosonic modes inside the super- 
fluid, while they still have negative energies that encour- 
age condensation. The mixing perturbation produces a 
new ground-state (with lower energy than the superfluid) 
as a quantum superposition of many different vortex ar- 
rangements and pushes away from it the other collective 
modes toward higher energies. The resulting low energy 
spectrum is qualitatively different than that of the trivial 
quantum Hall insulator because the ground-state is still a 
collective state of correlated bosons at some negative en- 
ergy where condensation is prevented only due to strong 
quantum fluctuations. 

We cannot predict the properties of the vortex liquid, 
and they are not universal. However, the perturbation 
theory requires that a vortex liquid be adjacent to any 
type-II superfluid phase of a two-dimensional quantum 
fermionic system at zero temperature with explicitly bro- 
ken time-reversal symmetry. How should this be inter- 
preted in the light of the fact that the naive perturba- 
tion theory does not correctly describe the lowest energy 
scales? Even though quantum fluctuations destroy the 
superfluid long-range order, they leave behind the bro- 
ken symmetries due to the vortex lattice. The vortex 
lattice potential energy (fT5)l is shaped by quasiparti- 
cle and Cooper pair mediated forces between vortices, 
which in turn are well defined at sufficiently short length- 
scales to be captured by the naive perturbation theory. 
Specifically, close to the pairing transition when the error 
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length-scale 8 becomes large in comparison with the sep- 
aration l m between vortices the naive estimate (|13[) can 
be trusted. Similarly, vortex mass is always dominated 
by the "high energy" degrees of freedom below the ultra- 
violet cut-off associated with the finite vortex core size, or 
the superconducting "gap" . Vortex is a local deformation 
of superfluid phases so that its dynamics can be reliably 
captured by the naive perturbation theory. Therefore, 
we can trust the estimate (|15|) as well and establish the 
existence of a vortex liquid at zero temperature. 

E. Renormalization group analysis 

The purpose of the following discussion is to elucidate 
the issues behind attempts to apply the simple theory 
([2]) to the analysis of realistic systems. A usual theoret- 
ical strategy is to justify a tractable field theory in the 
vicinity of second-order phase transitions where most op- 



erators allowed by symmetries do not significantly affect 
the macroscopically observable properties of the system. 
The fact that only a few operators are important to keep 
in calculations is what makes a field theory quantitatively 
useful. The model ([2]) does not have this useful property. 

We demonstrate using renormalization group (RG) 
that all arbitrary-range interactions between fermions are 
perturbations to @ that can shape phases and induce 
transitions. Physically, this comes from the fact that the 
bare fermion spectrum in (|2|) consists of dispersionless 
macroscopically degenerate Landau levels with energies 
e n = nu> c —fj/ (here we measure energy with respect to the 
redefined chemical potential fj/ = fj, — ui c /2). The lack of 
dispersion reduces dimensionality in RG and makes even 
arbitrarily weak perturbations extremely potent in lifting 
the degeneracy. We set up RG by generalizing ([2]) to d 
dimensions in the Landau level basis using the Landau 
gauge, and include additional allowed terms: 



S 



$t fr(o) $ 

2yj- «l:Px «1,«2 "2, 



J I " V ' n 1 n 2 J 

nm\rri2 

\- f dk xl dk x2 dq x , + f 1 

U2 2^ / o„ o_ o_ 1 mi...m 4 q x ) V mu k xl V m2 ,k as2 V m3 ,k x2 +q lc ' , P mi ,k xl -q x ( ^ • 



(16) 



We suppress spin a and flavor i indices for brevity. Since 
the bare fermion states are localized in the plane per- 
pendicular to the axis of rotation, it is appropriate to 
not rescale k x in RG. At zero temperature and density, 
fermion self-energy vanishes so that all renormalization 
comes from the boson field The RG flow equations can 
be calculated exactly to all orders of perturbation theory 
when U2 and the remaining omitted couplings are zero, 
since then the renormalization of $ involves summation 
of a geometric series of bare fermion bubble diagrams. 
Following the procedure in Refpi^, we find that under 
rescaling 

/l = e~ l r± t' = e- 2l T (17) 
ft = eW/'-i)^ $' = e W2-i)J$ 

the exact flow equations are: 
du' , du> r 

ir = 2 " if = ^ < 18 ' 

where a and b are cut-off dependent constants whose val- 
ues are not important for the present discussion. Appart 



from the Gaussian fixed point at // = u> c = v = g = 0, 
we can identify an interacting fixed point at fj, 1 = u) c = 
v = 0, g* = y/(3 - d/2)/(bN) which defines the unitar- 
ity limit in any dimension d, and shows how perturba- 
tion theory can be justified for large N. In addition to 
//, lu c and v, various couplings u n , which may include 
non-local short-range potentials and hence multi-particle 
collision terms, can be relevant at this fixed point. For 
n-particle scattering u n we find sufficient indication at 
the tree-level: 

-^ = [d+(2-d)n]u n + 0(u 2 ) (19) 

irrespective of the spatial dependence of interaction po- 
tentials in the plane perpendicular to rotation axis, since 
coordinates do not rescale in this plane. Therefore, all 
u n are relevant in d < 2, while in d > 2 they are relevant 
for n < d/(d-2). 

The same conclusions have been reached in the past 
using e = 6 — d expansions in the context of classical 
Landau-Ginzburg theories 29 . A formal advantage of the 
present approach is that the conclusions are justified in 
any number of dimensions as long as N is large, but also 
for not so large N as long as the u n couplings are small 
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and the density of particles is close to zero. 

The existence of infinitely many relevant directions 
in the parameter space near the Gaussian and unitar- 
ity fixed points implies the possibility that there could 
be many stable interacting fixed points that specify the 
properties of many different phases. The microscopic de- 
tails of interactions between particles can then very sen- 
sitively determine the concrete properties of a concrete 
system at low temperatures. We expect that the prop- 
erties of the vortex liquid in particular are not univer- 
sal. The vortex lattice melting transition lines are also 
not universal, being first order. The supcrfluid phase is 
generally amenable to phenomenological descriptions re- 
gardless of the microscopic details, but we can expect 
that the structure of vortex cores is system dependent. 

F. Conclusions 

We analyzed a simple model of neutral fcrmionic par- 
ticles in two dimensions with attractive interactions and 
explicitly broken time-reversal symmetry The model is 
routinely realized in ultra-cold atom experiments with 
fcrmionic atoms near a broad two-body Feshbach reso- 
nance, but it also captures some essential properties of 
much more complicated electronic systems. We assumed 
that the low-temperature superfluid is type-II so that it 
hosts an Abrikosov lattice of vortices at a fixed flux den- 
sity. Using a semi-classical perturbative expansion we 
found that a second order phase transition separates the 
superfluid from a normal phase of thermally excited un- 
paired fermions in the quantum Hall regime. We calcu- 
lated the critical temperature as a function of interaction 
strength and density and showed that in some circum- 
stances for weak interactions this phase transition can 
occur even at zero temperature. 

The unusual properties of the quantum Hall insulator 
at zero temperature were then exploited to argue that 
the second order superfluid-insulator transition obtained 
in the naive perturbation theory must be preempted by a 
first order transition due to quantum vortex lattice melt- 
ing. This conclusion followed from the analysis of quan- 
tum fluctuations related to vortex dynamics. We com- 
pared the potential energy of a vortex lattice with the 
kinetic energy of a vortex liquid and showed that the lat- 
ter always wins on the superfluid side close enough to the 
pairing transition. The vortex liquid is a strongly corre- 



lated normal phase which cannot be smoothly connected 
with the quantum Hall state of unpaired fermions. 

Finally, we argued that the properties of a vortex liq- 
uid in two dimensions at a finite flux density are not 
universal. Different materials can exhibit different nor- 
mal states resulting from vortex lattice quantum melting, 
including states with topological order such as quantum 
Hall states of Cooper pairs, or states with broken symme- 
tries involving density waves for example. At zero tempe- 
rature we expect that such states are insulators. At finite 
temperatures, below a crossover temperature, the normal 
states will be at least influenced by the proximate non- 
trivial zero-temperature phases, but we cannot rule out 
the existence of genuine finite temperature vortex liquids 
(depending on the broken symmetries of the zero tempe- 
rature insulator) which are separated from the disordered 
normal phase by a Kosterlitz-Thouless transition. 

Ultra-cold atom experiments have the potential to di- 
rectly explore the phenomena discussed in this paper. 
The target regime would be a combination of fast rota- 
tion and low particle density where a few lowest Landau 
levels would be occupied in the quantum Hall state at 
low temperatures. It is in this limit and on the BCS side 
of the Feshbach resonance where multiple isolated super- 
fluid phases could be found separated by normal states 
as the density of particles or rotation rate is varied. In 
order to faithfully represent our model the experiment 
would need to achieve this limit with a trap whose har- 
monic frequency is not much higher than the rotation 
rate. The phenomena predicted here might survive even 
in more confining traps, but the structure of vortex lat- 
tices and hence the phase boundaries would be affected 
by the trap in a manner that we cannot predict with a 
uniform model. Observing Cooper pairs at low tempera- 
tures in the normal states would be a first weak indicator 
of the proximity to a vortex liquid. 
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